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ABSTRACT 

Concise review of the basic properties of unitary matrix integrals. They are studied with the help of the three matrix models: the 
ordinary unitary model, Brezin-Gross-Witten model and the Harish-Charndra-Itzykson-Zuber model. Especial attention is paid to the 

i tricky sides of the story, from De Wit-t'Hooft anomaly in unitary integrals to the problem of correlators with Itzykson-Zuber measure. 
Of technical tools emphasized is the method of character expansions. The subject of unitary integrals remains highly under-investigated 

i and a lot of new results are expected in this field when it attracts sufficient attention. 

<N ■ 

2 ! 1 Introduction 

qq , Integrals over non-Hermitian matrices attracted far less attention during entire history of matrix model theory [1]- 
t-H ' [27]. This is unjust, both because they are no less important in applications than Hermitian models and because 
they can be effectively studied by the same methods and are naturally included into the unifying M-theory of 
matrix models [16]. Still, in the recent reincarnation of the matrix model theory [6]-[18] non-Hermitian models 
are once again under-investigated, and neither their non-trivial phase structure a la [6, 11], nor string-field-theory- 
like reformulation a la [12], nor combinatorial solution in the Gaussian (single-cut) phase a la [18] are explicitly 
analyzed - despite all these subjects are clearly within reach of the newly-developed approaches. 
i ■ The present text is also too short to address these important issues. Instead it is concentrated on the complement 

of these universal subjects: on methods and results specific to one particular - and most important of all non- 
t-H ■ Hermitian models - to that of Gaussian unitary ensembles. Importance of unitary integrals is obvious from the 
Yang-Mills-theory perspective: they describe angular - color - degrees of freedom, the complement of diagonal - 
colorless or hadronic - components, which are adequately captured by the eigenvalue Hermitian models [4] and 
jy-j \ by their GKM (Generalized Kontsevich Model [4, 21, 22]) "duals". As already mentioned, this does not mean 
(f) ■ that various kinds of eigenvalue techniques - from orthogonal polynomials [1] to Virasoro constraints [23] and their 
iterative [24, 11] or QFT-style [12] solutions - can not be applied to unitary integrals - they can and they are 
pretty effective, as usual [25, 26]. This means that of the main interest are questions of another type, concerning 
Q\ \ the correlators of non-diagonal matrix elements rather than those of traces or determinants. 

A powerful - and partly adequate - technique for the study of these aspects of unitary integrals is that of 
character expansions, and one of the main intrigues in it is the interplay between the two relevant groups: linear 
(GL) and symmetric (S). This interplay is now attracting a new attention because of the close link found between 
Yang-Mills theory - best represented for these purposes by lattice Kazakov-Migdal model [27]- [33] - and Hur- 
witz theory of ramified coverings of Riemann surfaces [34]- [40] - represented by the Hurwitz-Kontsevich partition 
function [35]- [38] and associated peculiar matrix model [36, 37]. 

This paper is no more than a brief introduction into these aspects of unitary integrals. Our presentation is 
mostly about the three basic unitary matrix models: 

• original unitary model [25], describing the correlators of traces of unitary matrices by the standard matrix- 
model methods, thus we touch it only briefly; 

• BGW (Brezin-Gross-Witten [41]) model, describing the correlators of arbitrary matrix elements, which was 
made - perhaps, surprisingly at that time - a part of the GKM theory in [32, 33] and was finally incorporated into 
the M-theory of matrix models [16] in [42]; 

• IZ (Itzykson-Zuber [43, 44]) model, a simplified single-plaquette version of the Kazakov-Migdal theory [27], 
describing arbitrary correlators with non-trivial but exactly-solvable weight (IZ action), - an old and difficult 
subject [29]-[31], which attracted some new attention recently in [14]. 

A lot of new - and much more profound - results are expected in the study of unitary integrals in foreseeable 
future, which - unlike that of Hermitian integrals - is in a rather early stage of development. This makes the 
presentation of this paper essentially incomplete and temporal - it is no more than a preliminary introduction. 
Still, we tried to concentrate on facts and formulas, which have good chances to remain in the core of the subject 
when it becomes more complete and self-contained. Instead some currently important - but temporal - results, 
i.e. those which will supposedly be significantly improved and represented in more adequate terms, are discussed 
in less detail and even unjustly ignored because of the space limitations. 
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2 Unitary integrals and BGW model 



2.1 The measure 

Unitary integrals are those over unitary N x N matrices U, UU^ — WU — I with invariant Haar measure [dU]. 
In variance with the case of Hcrmitian matrices H = H\ where dH = Ai^dHij, Haar measure [dU] is non-trivial 
and non-linear in U. If unitary matrices are expressed through Hermitian, one can express [dU] through dH, for 
example, 

U = e lH =► \ \SU\\ 2 = tr (U^SUU^SU) = tr (SH) 2 =► [dU] = dH, 

U= I - T ^^W\\ 2 = 4tr ( rJL= 6H-±- SH) m [dU] ~ (1) 

I-iH 1 11 \I + iH I-iH ) L J det N (l+H 2 ) V ' 

and so on. Alternative representation is through generic complex matrix M [32, 33, 14]: 

[dU] = <5(MM t - I)d 2 M = d 2 M j dS cxp (i tr S(MM^ - /)) (2) 

where S is auxiliary Hcrmitian matrix, used as Lagrange multiplier to impose the unitarity constraint on M. Below 
we often absorb i into S, what makes it anti-Hcrmitian, i.e. changes the contour of ^-integration. More than that, 
sometime one can absorb a non-trivial matrix factor into S, what completely breaks its hermiticity, still for the 
integral this is often no more than an innocent change of integration contour. 

In what follows we normalize the measure [dU] so that J[dU] = 1, i.e. divide [dU] by the volume Vn of the 
unitary group. This volume can be calculated, for example, by orthogonal polynomial method, see s.3.9 in the first 
paper of rcf.[4], and it is somewhat sophisticated: 

We also assume that the factor V^ 1 is included into the definition of Hermitian matrix integrals. 

2.2 Elementary unitary correlators in fundamental representation 

With this normalization agreement and the symmetry properties one immediately gets: 

ViAi) = J UiAiW] = jf5u6 jk (4) 

N 2 [SiiSiii'SjkSj'k' + Sii'Si'iSjk'Sj'k) - NiSu'Si'iSjkSj'k' + bii&vv&jk'bj'k) 
UijU^UlulAdU] = —± JV2 ( JV2_i) ( 5 ) 

It is clear from this example that the correlators with more entries become increasingly complicated and one needs 
a systematic approach to evaluate them. This approach is provided by the theory of BGW model [33]. 

Before proceeding to this subject in s.2.4 below, we note that these correlators are significantly simplified if we 
contract indices, i.e. consider the correlators of traces instead of particular matrix elements. Indeed, from (4) and 
(5) we get 



/ 
/ 



trtftrtft [dU] = 1, 
trC^tr^) 2 ^ = 2-S N1 , 
J tvU 2 (tTU^) 2 [dU]=S N1 , 
{tTU) 2 (tTU^) 2 [dU] = 2-S N1 (6) 



and so on (the case of N — 1 in these formulas requires additional consideration, see below). Such colorless 
correlators are studied by the ordinary unitary-matrix model [26]. 

Note that unitary correlators in (5) possess a non-trivial TV-dependence, with spurious poles at some integer 
values of N. This phenomenon is known as De Wit-t'Hooft anomaly [45]. Note also that in correlators of colorless 
quantities (6) these poles turn into peculiar non-singular corrections at the same values of N. 
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2.3 Prom the theory of unitary model [26] 

Partition function of the unitary matrix model is the generating function of correlators of colorless trace operators: 

Z v (t) = Jexp ( jr t k trU k ) [DU] (7) 

\k— — oo / 

It was also considered in [25]. Note that W = U^ 1 and conjugate [/-matrices are represented here by negative 
powers of U. Virasoro constraints have the form [26] 




d 2 d 2 



dt k dt n - k dt-kdtk-n, 

l<k<n x 

d 2 d 2 



Kk<n 



dtfrdtn-k dt-kdtk- 



Zu(t) = 0, n > 1, 



Zu(t) = 0, n > (8) 



The use of these Virasoro relations is a little unusual: there is no shift of any time variable and formal series 
solution is directly in terms of i's, not of their ratio with some background values of time- variables, which label 
the choice of phase, say, in Hermitian model [9]. The analogue of the Lq constraint states that 

X><lf = » < 9 » 

k— — oo 

and requests that monomial items in Zjj should have the form Y[k t^ k with X)/c>i ^ n k — Sfe>i ^ n -fc, i-e. allowed 
are tkt-k, tkh^-k-i and so on, but not tk or tkt-k-i with I ^ 0. With this selection rule, the first few terms of Z\j 
series are 

Z v = e Nt0 (l + a^ht-t + a 2]2 t 2 t- 2 + a lll2 (t 2 t_ 2 + t\t 2 ) + a^ntjt 2 ^ + . . . ) (10) 

They describe the correlators (6), and are dictated by the lowest Virasoro constraints with n = ±1,±2, but in 
somewhat non-trivial manner: 

•--""(j^s^+s&r)*- 

= t_i(-JV + JVai|i) +tii- 2 (-2a 1 | 1 +a 2 | 2 + 27Va 2 | 11 ) + t x t\ ( - JVai|i + a 2 |u + 2JVan|n) +... (11) 
Similarly, 

°--~''|£^ + ^ + |)^ = 

= i- 2 (-27V + iVQ !2 | 2 + 2a 2 | 11 ) +* 2 _ 1 (-ai|i+iVa 2 | 11 + 2a 11 | 11 ) + ... (12) 
Virasoro constraints require that each bracket at the r.h.s. vanishes, and this implies that 

«i|i = 1, 

(iV 2 - l)a 2 | 2 = 2(N 2 - ai,i) - 2{N 2 - 1) => a 2 , 2 = 2 - aS N>1 , 
(N 2 - l)a 2 |n =0 => a 2 | n = sS NA 

2{N 2 - l)an|n = {N 2 - lhn = N 2 - 1 a u \ u = \{2 - s5 N>1 ) 

(13) 

Thus Virasoro constraints do not define coefficients a 2 | 2 , a 2 \\\ and an|n unambiguously - at N = 1. In fact 
there is a single ambiguous parameter s in these three coefficients, and it can be fixed by direct evaluation of the 
correlators (6) at N = 1. Similarly, coefficients a with higher gradation level I = ^2 k>0 knk > 2 are not fully 
defined by Virasoro constraints at N = 1, . . . , I. 
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2.4 From the theory of BGW model [33] 

Partition function of the BGW matrix model [41] is the generating function of all unitary correlators: 
Z BGW = y>( jtc/+JC/t )[dt/] - j Je iS ^ MMl -^e trjlM+JM 'd 2 MdS = 

(2) f dS tl (jjt/ S )-trS S^l/S f dS tr (jjt S) _ trl/ s . . 

" J detS N J detS N 1 ' 

Obviously, Zbgw depends only on the eigenvalues of the matrix Jj\ i.e. on the time- variables tk = \tv ( JJ t )' £ . 

The model has two essentially different phases, named Kontsevich and character phases in [33], here we consider 
only the latter one, when Zbgw is a formal series in positive powers of tk- Virasoro constraints in this phase can 
be derived in the usual way: as Ward identities a la [23], but in this particular case they can be alternatively 
considered as direct corollary of the identity UW — I for unitary matrices, expressed in terms of i-variables. This 
identity implies that 

I j Zbgw = ~^^{J J^T^ 1 L n Z B GW - IZbgw (15) 

' 71=1 



0Jt dJ 



i.e. 

/ 

L n Z BGW = 

V 



E ktk 7f f + E T^T I Z BGw(t) = 6 nA , n > 1 (16) 



dtk+n fr* dt a dt b 

k=l a + b = n 

a,b>0 



are the usual discrete- Virasoro constraints of [19], familiar from the theory of Hermitian model, only with n > 1 
instead of n > —1 and with a modified r.h.s. at n = 1. 

Substituting a formal-series anzatz for Zbgw(^) into these Virasoro constraints one can iteratively reconstruct 
any particular coefficient, see s.3.2.1 of [33]: 

M>1 V fcl >...>fc„0 + + k M )\ 



and 



k 1 + ...+k M — 1 
CAr{fc Q } = C N {k a } Yl N 2 _ J2 

(18) 

1=0 



cjv(1) = N, 

c N (2) = -N, c N (l,l) = N 2 , 

c N (3) = 4N, cjv(2,1) = -3iV 2 , c N (l, 1, 1) = N(N 2 - 2), 

cjy(4) = -30AT, cjv(3, 1) = 8(2iV 2 - 1), 6^(2,2)^3(^2 + 6), 
cjv(2, 1, 1) - -6iV(iV 2 - 4), c w (l, 1, 1, 1) = A^ 4 - 8A^ 2 + 6, 



(19) 



The first two lines of this list reproduce eqs.(4) and (5). 
2.5 De Wit-t'Hooft anomaly [45, 33] 

Remarkable feature of (18) is the occurrence of poles at integer values of the matrix size TV. It deserves noting 
that this is the property of the character phase: nothing like this happens in Kontsevich phase of BGW model - as 
usual, if properly defined, Kontsevich partition function is actually independent of N [21, 33]. As we saw in s.2.3, 
these poles are associated with ambiguity of solution to Virasoro equations in the ordinary unitary matrix model. 
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In fact, DWH poles do not show up in correlators: numerators vanish at the same time as denominators, and 
singularity can be resolved by l'Hopitale rule. However, some care is needed. For example, in (5) one can consider 



J ul^ulfW] = 



2(N 2 -N) _ 2 N= i 



N 2 (N 2 -l) N(N+l) 

in accordance with (6). However, if one first take traces in (5), and then put N = 1, then one does not obtain (6): 

Anomaly of this kind is familiar from experience with dimensional regularizations, where taking traces and putting 
dimension equal to particular value also are non-commuting operations. 



2.6 BGW model in the M-theory of matrix models [42] 

As explained in [16] (see also [46] and [9] for preliminary results), partition function of every matrix model in every 
phase can be decomposed into elementary constituents, belonging to Kontsevich family of [21]. The most important 
building block is the original cubic Kontsevich tau- function tk, while BGW tau- function, made from Zbgw in 
Kontsevich phase, appears to be the next-important one (for example, it arises, along with tk, in decomposition 
of the complex matrix model. Both tk and tbgw satisfy the same continuous- Virasoro constraints of [23], only 
with n > — 1 in the case of tk and with n > in the case of tbgw- Another, compensating, difference is that 
instead of the second time- variable in the case of tk, for tbgw shifted is the first time: this allows the smaller set 
of constraints to unambiguously defined the partition function. See [42] for many more details of this important 
construction. 



2.7 Other correlators and representations 

A slight modification of BGW model is to make the power of det S in (14) an additional free variable. In this way 
we get the celebrated Leutwyler-Smilga integral [47], 

Z BG w{v\t) = J '(det Uye tT ^ u+Jut ^[dU] (20) 

which plays important role in applications and is widely discussed in the literature. 

Another modification involves the change of representation: in (14) partition function is defined as generating 
the integrals over unitary matrices in the fundamental representation of dimension N of GL(N), but one can 
instead consider an arbitrary representation R of dimension Dr(N). All correlators will be, of course, related 
to those in the fundamental representation, but in a somewhat non-trivial way, and they will be different. For 
example, instead of (4) in generic representation we have 

J U ab ul d [dU] = ^ (21) 

Moreover, different representations are orthogonal in the sense that 

8rr< (22) 

An important way to handle this kind of generalization is through the comprehensive theory of character expansions. 
The origins of this theory are in the fundamental works of G.Frobcnius and it is widely represented in modern 
theory, from pure mathematics to applied physics. As to references, the most we can do is to cite just a few basic 
textbooks [48] and just some related papers in the matrix-model literature [2, 3]. Our actual presentation will 
mostly follow [40], while applications of character theory to BGW and IZ models in s.4 will be discussed along the 
lines of [32, 33] and [49]. 



3 Basics of character calculus 
3.1 Three definitions of characters 

Character depends on representation R and on conjugation class of a group element, which can be represented 
through time- variables t\,t2, ■ ■ ■ or through auxiliary matrix (Miwa) variable A, accordingly we use two notations 
X R {t)=XR[X\. 
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Representations R of GL(oo) are labeled by Young diagrams, i.e. by ordered integer partitions R : \\ > X2 > 
. . . > 0, of \R\ - the number of boxes in the diagram, Xj = \R\, all Xj are integers. 

Time- variables form an infinite chain, which - if one likes - can be considered as eigenvalues of the infinite matrix 
from GL(oo). Miwa variables define an ./V-dimensional subspace in the space of time- variables, parameterized by 
the eigenvalues Xi of an N x N matrix X 1 in the following way: tk = jtr X k = j^2 { x k . Here and below tr 
without additional indices denotes a trace of an N x N dimensional matrix, i.e. a trace in the first fundamental 
representation of GL(N), associated with a single-box Young diagram R = 1 = □. The same matrix X can be 
converted into arbitrary representation R of SL(N), then we denote it by X^ . While X can be actually considered 
as independent of TV, i.e. defined as an element of GL (00), with N just the number of non- vanishing eigenvalues, its 
conversion into X^ actually depends on N , thus one should be a little more accurate with formulas, which involve 
X( R ) . However, for many purposes the TV-dependence enters only through dimensions Dr of representations R of 
SL(N), see eq.(30) below. In other words, while X^ and Dr essentially depend on TV, if formulas are written in 
terms of these objects, N often does not show up. Accordingly we suppress label N in X^ and Dr. 



There are three important definitions of characters: 

• the first Weyl determinant formula 

X*(t)=det(* Ai _ i+J -(t)) (23) 
expressing Xfl(*) through determinant of Shur polynomials, which are defined by 

exp(j>z fe ) =J2*k(t)z k (24) 

V k J k 

• the second Weyl determinant formula 

XR (t) = XR [X] = = d f* jX ' - (25) 



A(Y) det 



y x i 



where Van-der-Monde determinant A(X) = Yl^-ixi — %j) — dety x 1 ^ — det^ x\ 1 , 



i<j * 

• the trace formula 



XR {t) - X r[X] - Tr fl X^ (26) 



Actually, like a Janus, these characters have two faces. They are actually characters of two very different 
universal groups: GL(oo) and S(oo), where Sn is a group of permutations of N elements. Intimate relation 
between these two groups is already reflected in the definition of Miwa variables, expressing time-variables as 
symmetric functions of Y-eigenvalues. 

Important role in the character calculus is played by two additional functions: dimension of representations of 
SL(N) and Sn, associated with the Young diagram R: 

D R = dim R (GL(N)) = X R^tk = j^ =XR[I], N=\R\ (27) 

and 

d fl =-^dim fl (s(JV))=Xfl(tfc = 5fci), N=\R\ (28) 
Parameter d,R is given by the hook formula: 

dR = TT r4 r = dct 7^ " (29) 



and the ratio 



hook length i<i,j<|ii| (Xj + i - j)\ 

over all boxes of J 

Young diagram R 



Dr u (X i + N-i)l 

dR-U (N-i)l (30) 
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Note that d R , obtained by division of dimension by TV! = |i?|!, depends only on parameters Xi, while Dr contains 
an additional explicit dependence on N . 

One more important fact is that determinant of a matrix, i.e. a product of all its eigenvalues, is also equal to 
a character in the totally antisymmetric representation R = [1, 1, . . . , 1], 



ActX = X n...i[X], T&nk(X)=N (31) 

Moreover, integer powers of determinant are also characters: 

(detX) v = Xuu...A x l rankpf) = N (32) 

Note that when rank(X) = N the time-variables t k with k > N are algebraic functions of those with k < N, and, 
when restricted to this subspace, the characters with \R\ > N also become algebraic functions of characters with 
\R\<N. 

3.2 Sum rules and orthogonality relations 

Characters satisfy a number of "sum rules": certain sums over all representations (i.e. over all Young diagrams 
or over all ordered partitions) are equal to some distinguished quantity. So far the most important in applications 
were the following three sum rules: 

Sd«XflW=e tl =e tlX (33) 

R 



E XRiW) = cxp ( £ kt k t' k ) = Dct{mi _ x ^ x/) (34) 



vfe=l 



E J^Xr[X] X r[X'\ = J [dU] exp(trXC/X'C/t) (35 ) 

R 

Eq.(33) is a particular case of (34), because dn = Xfl(^fci)- Eq.(35) is actually a character expansion of Itzykson- 
Zuber integral, but it can be effectively used as a sum rule as well. 

A very important role in the formalism of character theory is played by orthogonality relation: 

XR(d)XR'(t) = S R , w for \R\ = \R'\ (36) 

Here one substitutes time-derivatives dk = = instead of t k in the arguments of Xfl(^) ~ a polynomial of 

its variables. For example, see s.3.3, 

d 



Mt) = t 2 + ^ =, x.(8) = -( W2 + w J, 



*? - r*s _ 1 ( d 92 



Xn(*) = -*2 + ^, =► Xii{d) = -[-— + 1 - s ) (37) 



and one can easily check that (36) works in these simple cases of \R\ = 1,2. It is also clear that orthogonality 
relation fails for \R\ ^ \R'\. 

3.3 Examples of characters 

For illustrative purposes we list here the first few characters and some other relevant characteristics of the lowest 
representations R of GL and S groups. 
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R 


AH 


rf_R = 

V nl m -i 1 


dim fl (S'| J j|) 

— 1 R|l r/ D 

— I J l-l - u iJ 


dim fl (SX(n)) 

AH \ fc / 


Vfl(2) = 

coeff t 2 t^ 1-2 
2dj? 


¥>i?(3) = 

coeff t 3 4 R '~ 3 
3d a 


¥>*(22) = 

coeff t 2 2 t[ R] ~ 4 
2 2 d R 




i 
i 


1 
1 


1 
1 


1 
1 


n 
u 


n 
u 


n 
u 


(1) 




1 


1 


n 











(2 s ) 






i 

2 


I 


n(n+l) 
2 




n 


n 


(1.1) 


+2 

- t2 + l 




1 

2 


1 


n(n — 1) 
2 


-1 








(3) 


f 3 

*3 + *2*1 + 

6 




1 

6 


1 


n(n+l)(n+2) 
6 


3 


2 





(2 ~\) 


*i 




1 

3 


2 


n(n 2 -l) 
3 


n 


_1 


n 


(1.1,1) 


f 3 

h - t 2 ti + -j- 

6 




1 

6 


1 


n(n-l)(n-2) 
6 


-3 


2 





(4) 


U + *3*1 + | + ^ + 




1 


1 


n(n+l)(n+2)(n+3) 


6 


8 


3 


24 


24 


24 




(3,1) 






1 


3 


(n— l)n(n+l)(n+2) 


2 





-1 




8 


8 


(2,2) 






1 

12 


2 


12 





-4 


3 








1 


Q 
o 


(n+l)n(n-l)(n-2) 




n 

u 


_i 




8 


8 


(1,1,1,1) 


t| t 2 t 2 
-t 4 + t 3 h + -j - + 




1 


1 


n(n-l)(n-2)(n-3) 


-6 


8 


3 


24 


24 


24 




(fc) 

\ / 






J_ 


1 

= |i2|! 


(n+fc-1)! 


fc(fc-i) 


fc(fc-l)(fc-2) 


fc(fc-l)(fc-2)(fc-3) 






fc! 


(n— l)!fc! 


2 


3 


8 
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3.4 Example: eq.(26) for GL(2) 

Construction of from a given X is a little tricky. Actually, in explicit form eq.(26) is a relation for unitary 
matrices: 

XR [e A ]=Tv R e AW (38) 
and relation is simple between A and A^: the n x n and D R x D R matrices, both expressed through the same 



set of parameters u a : 



n"-l 



A = i"£u a T^ and A™ = i ]T u a T™ 

a=0 a=0 

Here are generators of the GL{n) algebra in representation R, and R = 1 = □ is the first fundamental 

representation. 

For example, in the case of GL(2) 

a + d b 

c a — d 



A = 



and we have: 



Xi [e A ] = ti = Tr Dl=2 e- 4 = Tr 2 jexp 



a + d b 

c a — d 



X2 [e A ] - t 2 + ^ - -Tr 2 e 2A + - (Tr 2 e A ) 2 = Tr D2=3 ^ exp 2c 2a 2b 

1 ~ 1 < c 2a-2d 

Xii [e A ] = ~h + | = -^Tr 2 c 2 ^ + I (Tr 2 e ^) 2 = e 2 «, D n = 1 



1 

~ 3 


T¥ 2e ^ + 1(TV 2 


e 2A ) (Tr 2 e A ) 


+ J(^) 8 








/ 3a + 3d 


b 


\ 




=4 < 




3c 


3a + d 2b 







exp 





2c 3a — 


d 36 


> 






V o 


c 


3a — 3d J 


> 



X21 [e A ] = -t 3 + | = -^Tr 2 e* A + 1 (Tr 2 e^) 3 = Tr D21=2 jexp ( 



3a + d 6 
c 3a — d 



Xin [e A ] = t 3 - t 2tl + | = i T r 2 e 3 - 4 - 1 (Tr 2 e 2 ^) (Tr 2 e A ) + \ 



Tr 2 e A ) 3 = 0, D m =0 



Note that: 

• matrix A does not actually need to be anti-Hcrmitian, above relations are actually true for arbitrary A; 

• normalization of the central U(l) generator includes a representation-dependent factor \R\; 

• in general, for representation of the size rxrof GL{2) we have: 



TV = 



/ 1 

2 





\ 

/ r - 1 




To = 





V o 





r-3 






o \ 




r-2 



T_ = 



r- 1 
j 






\ 




3-r 
1-r / 



with commutation relations [T+,T_] = T 3 , [T 3 ,T±] = ±2T±. 



/ o 

r - 1 
r 



v o 

To = \R\ 














\ 
















- 2 


















2 
















1 


o / 


/ 1 










\ 





1 


















1 





V o 










1 J 



(39) 
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3.5 Characters and symmetric functions [32, 49] 

Many application of characters are based on the fact that they are symmetric functions (of eigenvalues Xi of a 
matrix X), and simultaneously arc cigcnfunctions and eigenvalues of certain distinguished operators. 
In the space of symmetric functions there are three "obvious" sets of homogeneous generators: 

time - variables p k [X] = = irX k , k > 0, p k [X] = kt k [X], (40) 

i 

Shur polynomials Sk[X] = ^ x^ . . .Xi k , k > 1, s [X] = 1 (41) 

»1<...<H 

and 

elementary symmetric polynomials <Jk[X] = ^ x il . . .x ik , k> 1, <Jo[X) = 1. (42) 

ii <...<»*, 

The corresponding three generating functions are 

OO j 

P( z )=Y / P k [X}z k = trj—^ (43) 

fc=0 



oo - 

*(*) = = n ( i - - ^T^y W 

and 

oo 

£(*) = ^ = II^ 1 + = dct ( 7 + zX ) (45) 

fe=0 i 

These three generating functions are related as follows: 

Siz)=e^(±t k [X\A=^ ( ^ (46) 



where [X] = ktk [X] , so that 



P(z) = tTl + z-^-logE(z) (47) 



and 

£(z) = (48) 

The sets {tfc}, {st} and {er/c} are generators in the space of symmetric functions: any such function is a poly- 
linear combination of generators. Linear bases are formed by products of these generators and are labeled by 
Young diagrams R. An obvious choice of basis is 

3 3 3 3 \ 3 J 

In this way one can define the "i?-power" of any semi-infinite sequence {pi}, not obligatory from the time- variables. 
Sometime we need the same basis with a different normalization: 

R 

p^S = ?— with z R = TT \j m i (50) 
zr ■ lx 

3 

Note that widetilde here is not exactly the same as tilde in orthogonality relation (36). 
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A less trivial linear basis is provided by characters, made from the generators {s k } with the help of (23), 

Xr = det (s x j-j+i) (51) 
y 

Motivation for this choice comes from the following expansion rule [49]: 

II ( £ J = SxflW • det (A Xj -j+i) (52) 
which is related through the second Weyl formula (25) to the basic determinant expansion [32, 49] 

(\ N N 

^A tk B 3k = {-) p YY}ZA lk B P(l)k = J2^ p H Ai ^ Bp ^ = 
k ) PeS N i=l k fci,fe 2 ,...,feiv P »=1 

= E (E(-) P (-) P 'II^'Wfci B ^Wfci ) = E detA ife .detS ife . (53) 

fei>fe 2 >... \p,p' i=i y fc 1 >fe 2 >... y y 

which generalizes det(AB) — det A ■ det B from square to rectangular matrices. Indeed, with kj = Xj — j (note 
that from Ai > A2 > . . . it follows that k\ = \± — 1 > k?. = X2 — 2 > . . .) and with (25) substituted for xr[X], the 
r.h.s. of (52) has exactly the same form as the r.h.s. of (53). Therefore 

_ ,( 25 ) det ij ( X i J ~ 3 ) det iJ ( A ^-j+i) 

YxR[X}-det(A Xj - j+i ) W J2 dc 7 t x -i = 

r j Ai>a 2 >... aez v x J 

(53) dctjj (gfc A fc+i ^) = dct tJ - (A{ Xj )x^) = j-j-^-j (M) 

£„• 

with A{x) = Y^k=o A k x k and we deduce (52) from (53). Note that summation limits are automatically adjusted in 
appropriate way. 

3.6 Characters as eigenfunctions and eigenvalues [36, 38] 

Characters Xi?[^] are the common eigenfunctions of the infinite set of commuting operators W(R) [38]. These 
operators are in one-to-one correspondence with the Young diagrams (ordered integer partitions) and are made 
from the GL generator D t j = X ik d ^ - in two steps: first define a semi- infinite sequence of commuting operators 

D k = : tr D k : and then use the general rule (50) to introduce 

W(R) = :m: = :]J-^(tvD k r: (55) 

k 

Additional ingredient is the normal ordering, which implies that all X stand to the left of all d/dX. According to 
[38] these operators form the center of universal enveloping algebra, 

W(R),W(R') =0 VR,R' (56) 

and characters are their eigenfunctions 

W{R) X r> [X] = <p R , {R)xw [X] WR, R' (57) 

The corresponding eigenvalues <pr{A) are in fact characters of symmetric group S, and they can be defined as the 
coefficients of XR expansion in p- variables, 

XR(P)= E d RVR{^)p A (58) 

|A| = |fl| 

The sum here is over Young diagrams A of the same size (number of boxes) as R. For A > \R\ these eigenvalues 
vanish, while for |A| < \R\ they are defined as 

k \R\-\A\ 
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i.e. if A does not contain 1 entries at all. In particular, for such A we have ipr(A) = <p R (A, 1, . . . , 1). In other 

|AHA| 

words, if A contains exactly k unit entries, then 

|fl|-|A| 

Important in applications is also inverse of (58): 

P*= d RX R(p)cpB.(A) (61) 

|fl| = |A| 

Note that p A appears at the l.h.s. instead of p A at the r.h.s. of (58). Note also that the both formulas are written 
in terms of p rather than t- variables, pt — ktk- 

Operators W{R) and thus their eigenvalues form an interesting commutative associative algebra [38, 40] 

W(R 1 )W(R 2 ) = '£C R l3 lR W( R z) => <Pr(Ri)<Pr{R*) = £ C rIr^r{Rz) Vi? (62) 

4 The theory of IZ integral 

The very important application of character theory is the theory of Harish-Chandra-Itzykson-Zuber integral [43, 
44]. As mentioned in the Introduction, it is the simplified version - or the main basic block - of the Kazakov- 
Migdal model [27]- [33] of lattice gluodynamics and, actually, of many other matrix- model-based approaches to 
non-perturbative dynamics of gauge theories with propagating (particle-like) degrees of freedom. It is also a non- 
trivial example of group integral, which satisfies the requirements of Duistermaat-Heckman theorem [50] and can 
be treated exactly by the methods, now known as localization technique (with Nekrasov's instanton calculus [55] 
in Seiberg-Witten theory [56] as one of its remarkable achievements) . 

HCIZ integrals can be handled by a variety of methods, starting from the Ward identities of [44]. We present 
here the character-based approach [32, 49] and a closely related GKM-based technique of [29]- [33], recently- revived 
in [14] and applied to the old (and still unsolved) problem of arbitrary correlators with the HCIZ weight [30, 31]. 

4.1 Itzykson-Zuber integral and duality formula 

IZ integral is a unitary integral with non-trivial weight of the special form: 

J IZ (X,Y) = J e trXUYU '[dU] (63) 
It depends on the eigenvalues of two N x N matrices X and Y through the celebrated determinant formula [44] 

JlZ ~ A(XjA(Y) = A(X)A(Y) 2. H ex P (^W) ' ( 64 ) 

where A(X) — Y\ i< j(xi — Xj) — det^ x^ 1 and A(Y) are the two Van-der-Monde determinants. 

The standard way to derive (64) is to observe that Jjz depends only on the eigenvalues of X and Y and 
substitute a function of eigenvalues into the obvious Ward identities: 

{ tT {dX^) -^ Yk ^ J iz(X 1 Y)=Q 1 jtr(^r) - trX fc | J IZ (X,Y) = 0, k > (65) 

what reduces the problem to Calogero-Dunkl-like equations, since, for example, 

( d V d 2 1 f d d \ 

^{dX^ =1^0^ + 1^^— ^X^-dx-J (66) 

when it acts on a function of eigenvalues. One can then check that (64) is indeed a solution. This method is 
rather tedious and the answer is in no way obvious, at least for people unfamiliar with Calogero-Dunkl equations 
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and associated theory of "zonal spherical functions" [57]. Instead it has non-trivial generalizations, a recent one 
[58] is to the theory of integral discriminants: the homogeneous non-Gaussian integrals - an important branch of 
non-linear algebra [59]. 

Direct way from (63) to the answer (64) is provided by Duistermaat-Heckman approach [50], which claims that 
an integral invariant under an action of a compact group and with appropriate relation between the classical action 
and quantum measure - conditions, which arc trivially satisfied by the group-theoretical integral (63), - is exactly 
given by its full quasiclassical approximation: a sum over all extrema of the classical action: 



/ 



O : 5S(0o)=O ^ ' 

As explained in s.3.2 of [5], see also [14], the equations of motion for the action tr XUYU\ 

[X, UYtf] = 0, (68) 

after a unitary transformation which makes X and Y diagonal have any permutation matrix, U = P, as their 
solution (more solutions arise when some eigenvalues of X or Y coincide, but this does not affect the answer). 
Then Duistermaat-Heckman theorem (67) implies that 



J IZ (X, Y)~Y1 A(X) A(Y) CXP ( ? XiVp W ) (69) 



where Van-der-Monde determinants arise from the pre-exponential factor and the sign factor (— ) p naturally 
arises from the square root and distinguishes the contributions of minima and maxima to the quasiclassical sum. 
Dusitermaat-Heckman method, known also as localization technique, got a variety of applications since it was first 
discovered in the study of HCIZ integrals. A new promising application should be to the 3-algebras and associated 
BLG theory [60] , where this is still the only available direct approach to defining a substitute of matrix models - 
since no direct counterpart of matrices is yet known for 3-algebras. 

Of certain importance is also the following representation of (63) as a Hermitian-matrix integral [29, 33, 14]. 
Immediately from the definition (63) we get: 

JlZ ( = fe^ XUYU \dU]^ ( d 2 MdSe trS(I-MM^+trXMYM^ = f ^ S dS 

J J J det N 2 xN 2(S®I -X®Y) v ' 



One can easily convert to the basis where, say, Y is diagonal, Y — diag{j/fe}. Then 

s trS 

t tet(S-y k X) J I\ k dct(S-y k I) J " l »^~> Y[ uk { Sl - y k ) 



where determinants and traces are now over N x N matrices and arrow transition results from the change of 
variables S — > SX (which implies a change of integration contour over S) . The last step is transition to eigenvalues 
of S, with angular integration given again by the same Itzykson-Zuber integral. According to [14] this duality 
relation provides an effective recursion in the size N of the unitary matrix. 



4.2 Character expansion of the HCIZ integral 

From character calculus we can immediately derive an alternative formula, called character expansion: making use 

(33) 

of (33) one can express e tr * with ^ = XUYW as a linear combination of characters, e tr * = daXR^) and 

after that substitute the average of character over the unitary group: 

J Xfl [*] [dU] ( = 6) J Tr R * R [dU] ( = 6) J Tr R (x R Uy R tf)[dU] ( = ] ±- Xr [X} X r[Y] (72) 

Thus 

Jrz = [e«*[dU] ^ J> f XR m [dU] ^ 2 d *X*W**M (?3) 

•* R R R 

This elementary but powerful formula can be found, for example, in [49] and - with erroneously omitted (i^-factor 
- in [3]. For its recent application to exact evaluation of correlators in Gaussian Hermitian model see the last 
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papers in ref.[17]. Similarly one obtains a character expansion in a more sophisticated case, with two terms in the 
exponent, known as (generalized) Berezin-Karpelevich integral [51, 49]: 



>BK 



[dU] ( [dV] cxp (XUYV* + VYU^X) = ]T d R d R , f f Xr {XUYV^) X r'{VYU ] X)[dU][dV] = 

•* R,R' 

( = E4/ [Tr R (x R Uy R Vl)Tr R (x R Vy R tf)[dU}[dV} ^ J2^Xr[X 2 ] X r[Y 2 ] (74) 

R •* J R R 

because, as a consequence of (22) we have, in addition to 

Jxr[U}xr'[U] [dU]=6 R , R , (75) 

For applications to two-matrix models see, for example, [3]. 

Representation (64) follows from (73) if one expresses characters at the r.h.s. through the eigenvalues with the 
help of the second Weyl formula (25): 

(73) d R XR[X]xR[Y] (25) 1 d R N-j+Xj , , N-j+Xj 

JlZ = § Dr- = AfW)^^ 1 - = 

\ i +N-i=k i n^i(-^-*) ! \ - dctjj- x\ ] det tJ - yl 3 _ I jj \ det aj e x w 

l^i U .kA 1 11 I A( X)A(V) 1 ' 



since, as a particular case of (53) 



]T ( J] Oj ■ ) det ^ det y*' = ^ E E (M^P'w) ^ = d 2 ct (77) 
fci>fe 2 >... / y lJ ' p,p' fc. 



for A(x) = J2k a kX k , see [49]. 

Similarly from (74) we obtain another celebrated formula [51, 49]: 



I i ; 2fc,- , . 2fe,- / AT-1 



(74) x ^ d|xp[^ 2 ]XPy 2 ] = (n,7=o J 1 ) x ^ detjj x- " 3 detjj (77) / TT •, ] dct i3 1 {2x i y j ) 



where Bessel function Zo(2x) = J2 k pj. This is a kind of a transcendental generalization of the IZ integral. 
4.3 Reductions to BGW model and Leutwyler-Smilga integral 

Similar character expansions exist, of course, in the simpler case of BGW model. However, since integral (14) 
contains two items in the exponent, it is in fact closer to the transcendental Berezin-Karpelevich integral, than to 
the IZ integral. Instead of (74) we get for (14): 

Zbgw = f e^^ u + Jur \dU] = £ d R d w j XR [J^U] Xw [JUn [dU] (2 i } ]T ^XR[JJ"\ (79) 

■> R,R' R R 

To make the next step [49] we now need not only (30) but also the hook formula (29): 

J 2 FT^ 1 7'' dot T 2k i 1 



(53) 



R 

'AT-1 



D * A[JJt] ^ ( n .,.,) TT ■ - JV + i)! 



det« (jf-ili-sVJij) (n-i \ det« (j/- 1 2 J -i(2J i )) 
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where Jf are eigenvalues of Hermitian matrix JJ^ and the Bessel function 

oo 2/c+s 

^) = E^)T=^^) (81) 

Leutwyler-Smilga integral (20) can be handled in the same way, but one should make use of appropriate 
generalization of (33) and (29): 

{d etUre*v= £ X^AWte 1 (82) 

A!>A 2 >...>0 J \ J J) 

From this formula we get [47, 49] 

Z BG w(v,t)= f(detUre^^ JU %U} = ^^ £ ^M.et - 1 — - = 

J A t JJt i fcl >fe>... (n^ii) y fc-^ -"+*)' 

N-i \ dety (j?'- 1 J v+J -_ 1 (2J i )) 
4.4 Correlators of \I/-variable in IZ integral [40] 

Character expansion allows one to evaluate arbitrary correlators of traces of * = XUYU^ and its powers in IZ 
integral. 

For A = |<5i > (5 2 > . . . | = | m, . . . , m, ... , 1, . . . , 1 j of the size | A| = J2k $k = J2 m ma m we have from (61): 
* A = (trvI/)^(trvI/ 2 )^...(tr* m ) a ™ ( => z A £ d flX fl (*)</>fl(A) (84) 

|fl| = |A| 

Instead of (73) we can now write 

/tV^^y^^'^^l £ ^ fl (A + n-l)/x«W[^] = 

^ n=0 n ' n |R| = |A|+n 

= £ m - w' ^ y "( A + (|fi| ' |A|) • < 85 » 



Note that sum over n is actually traded for the sum over all R, then only the contribution with n = \R\ — |A| is 
picked up from the sum over n. If A does not contain any unit entries, then Za+(|.r|-|A|)-i = (\R\ — |A|)! za and 
this sophisticated formula simplifies a little: 

f^e t ^[dU}=Y / ^^R(^ + (\R\-m^)xR[X] X R[Y} if HA (86) 

R R 

As a byproduct of this calculation we obtain a possible proof of (57) [40]. For this purpose we act with 
the operator W'(A) = D A on Jiz- Since operator L>ij = J2k Xik d xik converts = (XUYU^)ki into itself, 
Dij^kl = Sjk^iU we have: 

W(A)J IZ {X,Y) = W{A) J e tr *[rf[7] = J *^e tr * [dU] (87) 



and now it remains to substitute (72) for Jjz at the l.h.s. and (85) for the average at the r.h.s. Picking up the 
coefficient in front of in sums over R at both sides we obtain: 

W(A) XR (X) = *n?-V\fz y *( A + ilRl ~ |A|) " O**^ (6=0) f^ A )XR[X] (88) 
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The last equality is trivially true when A does not contain unit entries, 1 ^ A, since in this case 2a+(|.r|-|A|)-i = 
(|.R| — |A|)! za and i^^^A+(|_R| — |A|)-1^ = ipr(A). If A contains k unit entries, then ZA+(|ii|-|A|)-i = ^ R ^^ +k ^' z A , 

while ^(A + aai-IAD-l) (6 = 0) ^J^^(A). 

Coming back to (85), it provides a character expansion for generic correlator of ^-traces. Of course, this is 
only a half-stuff formula: an analogue of the transformation (76) is still needed to convert it into a more explicit 
expression. Unfortunately, in any case the correlators of ^-traces are not the most general ones needed in HCIZ 
theory, and evaluation of other correlators remains an unsolved, underestimated and under-investigated problem 
to which we devote a separate short section. 



4.5 Pair correlator in HCIZ theory 

Evaluation of unitary integrals with HCIZ measure became an important scientific problem since the appearance 
of Kazakov-Migdal model [27] of lattice gluodynamics and was originally addressed in [29] and [30, 31]. Given the 
Duistermaat-Heckman nature of HCIZ integral (69) [5, 14], this problem is also a natural chapter in development 
of "localization" techniques and BRST-like approaches [50]- [57]. Of course, the problem is also related to generic 
studies of non-Gaussian integrals in the theory of integral discriminants [58] , which is an important part of the newly 
emergent field of non-linear algebra [59]. Unfortunately, the subject of HCIZ correlators is still underestimated and 
remains poorly investigated - not a big surprise, given the lack of knowledge about the ordinary unitary integrals, 
overviewed in s.2 above. Still, the interest is returning to this field, in particular the relatively recent papers [14] 
provide a new confirmation and serious extension of the old result of [30] . This section provides nothing more than 
a sketchy review of the present state of affairs. 

Correlators can be in principle evaluated directly - by choosing one's favorite parametrization of unitary matrices 
and explicitly evaluating all the N 2 integrals. A clever and convenient parametrization [31] makes use of Gelfand- 
Zeitlin variables, which are unconstrained (free), it also arises naturally in recursive procedure in TV of [14], which 
starts from duality relation (71). Though seemingly natural, this approach does not provide any general formulas, 
at least in any straightforward way. As usual, a more practical approach is the "stringy" one [53] - that of the 
generating functions. 

The generating function of interest in HCIZ theory is 

Z IZ = J e trXUYUUtrjW+trJU^ dU] P) 

- J d 2 MdScxp (trS(I - MM^exp (tr XMYM^ + tr J f M + = 

/e tvS dS ( / \ _1 \ 

det^W* ® J - X S r) ^^ J \ S ®I-X®y) J) = 

= / U k ^(S-y k X) CXP (*» + J ^ S - X ri« J ») } ^ 

- as direct generalization of (70) and (71). Integral over M is Gaussian and Wick theorem is applicable, non-trivial 
part is evaluation of the S integral. In can be interesting to include also the Leutwyler-Smilga factor (det U) v into 
the integrand, but - as we saw in s.4.3 - this leads to non-trivial modifications. Unfortunately not much is known 
about this function even at v = 0. There are no serious obstacles for going beyond the first non-trivial correlator 
(UijUl^iz, still it remains the only well-studied example. 

This pair correlator simplifies considerably if we consider it in the basis where X and Y are both diagonal - 
then [29] 

J UijUl [dU] = M l0 5 k0 8 l i (90) 
According to [30], in this basis the generating function 

J2 "iMijbj = J2 l^fcxp r£x kVl \U k A [dU] = A(x) X A(y) Yl (-) P e x > v ^>. (91) 

— 1 i,j=l \ k,l J P(zSn 
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/ 1 



det 



JV-l 



n-1 



EH" E 



X 

V ah 



1 \ 



X in+1 
din+1 ) 



det 



/ 1 

yp(*i) 



y P(ii) 



1 \ 

yp(i„+i) 



y P(in+l) 



\ b P{ll) . . . &p( in+1 ) ) 



n=0 



l<ii<...<i n+ i<N 



( 1 



det 



1 \ 



det 



/ 1 

yp(h 



\ x i ••• < +1 / 



1 \ 

yp(i n +i) 



V ypin) ■■■ yp(i n+1 ) J 



( 



A(X)A(Y) 



E\ - {a% - aj)(bp(i) - bp(j)) 
aibp(i) - > 
i<^< N _ x o)(yp(i) - vpu)) 



i=l 



V 



E 



1 1 1 

Xi Xj X k 
di dj d k 



1 1 1 

yp(i) yp(j) yp{k) 

bp(i) 6p(j) &p(fc) 



\<i<j<k<N 



(xi - Xj)(xj - x k ){ Xl - x k )(y P (i) - yp(j)){yp(j) - yp(k)){yp(i) - yp(k)) 



+ . 



This sophisticated formula was transformed in [14] into a more elegant form: 

1 f, „ , /„ 1 „ 1 



A{X)A{Y) 



res jdet£-det 



■E: 



X-u-I Y-vI 



(92) 



v =Vj 



where matrix E is defined as Eij — e XiVj . According to this residue formula one should pick up the coefficients in 
front of ^(u — Xi)(v — yj)^j in the sum 

1 



A(X)A(Y) p £ ( )PI ? 



e XkVp( k ) l 



(u - x k ){v - y P{k) ) 



(93) 



over permutations. The first term det E in (92) does not contribute to the residue. There will be two different 
kinds of distributions: one for permutations P with the property j = P(i) and another - for all other permutations: 



Mi 



1 



A(X)A(Y) 



E Eh j 



pgSjv p 



N ( 1 

5PWj ]ji \ ( Xi - Xk ^ V 3 - yp ik)) 



+ 



i 1 &P{l) ' 3 ) {Xi-Xp (j) ){ yj 



yp(i)) 



N 

n 

k^LP-Hj) 



1 - 



(xi - x k )(yj - y P{k) ) 



N 



(94) 



k=l 



In particular, for N = 1 we have a single term of the first kind and Mn = e XlVl . 
For N = 2 the terms of both types contribute: 



M 



1 



li 



M 12 = 



x\iy\i 
1 



1 



£122/12 



^C\y\+x 2 y 2 _|_ 



1 



1 



„x 1 y 1 +x 2 'L 



x\2y\2 yx\iy\2 



1 



x\iy\i 
1 

x\iy\i 



x 1 y 2 +x 2 yi 



x 1 y 2 +x 2 yi 



M 2 2, 
= M21 



(95) 



what is exactly the same as (91), where just the first two terms survive at the r.h.s. for N = 2. Of course, 
Ma + M12 = Jiz- 
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In general, equivalence of (91) and (92) follows from the fact that for a.j = l/(u — Xi) the determinants in the 
numerator in (91) turns into 



/ 1 



det 



\ 



„"-l 



c in+l 



I X ... 1 \ 



det 



n+l 



/c=i 



U — Xi 



(96) 



- simply because the numerator should have dimension n(n — l)/2 and be antisymmetric in all x ik . Therefore (91) 
implies that 



n+l 



1 



I n=0 l<u<...<i„ + i<JV fe= 



L x («-*i|b)(«-I/P(i»)) 



AT 



V (-) p e XiVp ^ < 1 — TT 11 - 

A{X)A(Y) p ^ > \ l\\ (u- Xi ){v-y P(i) ) 



det £- dot (E-^E^ 
A(X)A(Y) 



(97) 



Note that taking a limit to (4) in this formula is a somewhat tricky exercise - this is not a big surprise, since 
even Jiz{X, Y = 0) = 1 is not a fully obvious corollary of eq.(64). Much simpler is to check that (97) is consistent 
with V. \Uij\ 2 = 1, i.e. that V. M%j = Jiz for any i. Indeed, 



^ Ma =^ f dv 



det 



A(X)A(Y) 



(98) 



- / ^ (l - det (l - ^-E^-eA )=J Iz f tr ^-E^-E^ = J IZ tr -J— 

A{X)A(Y) \ \ u-X v-Y J J u-X v-Y u-X 

A character decomposition of pair correlator can be easily obtained similarly to eq.(72), provided one knows 
appropriate generalization of (21) and (22), namely 



J u R u R u R ,u RI [du] 



(99) 



- a generalization of (5) to arbitrary representations. Actually, the case when R' is the fundamental representation 
would be enough for decomposition of (91). 

Like HCIZ integral itself, the pair correlator satisfies a kind of duality relation - a direct generalization of (71): 
from 



/ UijUlidU] 



(89) 



Ojk 



1 



z tTS dS 



S-Xy k J u ILdet(S - Xy k ) 
it follows, after the change of integration matrix-variable S — ► XS and after diagonalization of X and S that 

M t3 (X,S) A 2 (5) U p ds P 



(100) 



Mij(X,Y) 



{Si-Vj)Xj Ylp,q(Sp ~Vq) 



(101) 



Similar relations can be straightforwardly deduced from (89) for higher correlators in HCIZ theory. 
For generalizations to other simple Lie groups see, for example, [7] and [14]. 



4.6 U{2) example 

To illustrate the use of duality relations we demonstrate explicitly how (101) is satisfied by (91) and (92) in the 
simplest non-trivial case of N = 2. 

We begin with an even simpler story: duality relation (71) at N — 1. The statement is simply that 



e xs ds 
s-y 



(102) 
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and it is clear that integral should be understood as taken around the pole at s = y. 
For N = 2 the same (71) becomes slightly more involved: 

e x iyi +x 2 y 2 _ e x iy2 +X2Vi (g9) f f e x 1S1 +X2S 2 _ e xi S2 +X2Sl sl 2 dsids 2 



£122/12 J J x 12 s 12 (si - vi){si - y 2 )(s 2 - yi){s 2 - 2/2) 

! ' f ( e xis 1+ x 2 s 2 _ e x lS2 +x 2Sl j f 1 1 \ f 1 1 \ Sl2dsids2 = 



X\2Vl2 J J " ' \Sl~yi S1-JJ2J \S2-yi S 2 -V2 



1 / x p xiyi+x 2 y 2 _ p x 1 y 2 +x 2 y 1 

—o-[ - yi2(e XlVl+X2V2 - e Xiy2+X2V1 )- y 2 i ( e Xiy2+X2yi - e XlVl+X2V2 ) ) = -2 (103) 

Xl2y( 2 V 1 £122/12 



Note that from the four possible poles, (si,s 2 ) = (2/1,2/1), (2/1,2/2), (2/2,2/1): (2/2,2/2) the first and the last one do 
not contribute because of the S12 = Si — s 2 factor in the numerator of the integrand. 
Coming back to (101), in the case of TV = 1 it states simply that 



e xs ds 1 de x 



x(s — y) 2 x ds 



= e xy (104) 

s=y 



For N — 2, the calculation is much longer. We present it only for one particular (out of four) non-vanishing 
correlator (95), 

Mu(X,Y) = 2^2- ((1 - x 12 yi2)e XlVl+X2y2 - e Xiy2+X2y A (105) 

With the help of the same tricks as in (103) we can write the r.h.s. of (101) as 

s\ 2 ds x ds 2 (I I \( I 1 \ \ (I s 12 \ _ X1S1+X2S2 + J_ oX1S2+X2S1 



(si -yi)xiyf 2 \S! -yi s 1 -y 2 J\s 2 -yi S2-2/2/I \x\ 2 x 12 J x 12 
dsids 2 { -7 1 rrr + 



( X 12 


«12 




Xl2 


1 


*)} 


si - 





x\x\ 2 y{ 2 J J | (gi - yi) 2 2/i2V s i-2/i «i - 2/2 / J V s 2 - 2/1 s 2 - 2/2 

•((1 - s 12 x 12 )e XlSl+X2S2 - e XlS2+X2Sl ^j 

This time all the four poles can contribute and cancellations for two unwanted ones are a little less trivial: for 
(si,s 2 ) = (2/1,2/1) we have 

(-{ Xl - x 12 ) +x 2 + — (1 - 1)) e^+ X2 ^ = 0, 
V y 12 I 

while for (s 1 ,s 2 ) = (2/2,2/2) 

— (1 - 1) e (xi+X2)v2 = 
2/12 

To simplify analysis of these formulas we underline the terms, which come from the double pole at s\ = y\ and 
produce derivatives of the integrand over s\. Non- vanishing are contributions from the two other poles, from 
(si,s 2 ) = (2/2,2/1): 

— — ( (1 + £122/12) e xlV2+X2V1 - e xiyi+X2V2 ) 
2/12 v > 

and, finally, from {s 1 ,s 2 ) = (2/1,2/2): 

(xi(l -xvtyii) -X12) e XlVl+X2V2 - x 2 e XlV2+X2Vl - -j- ((1 - x 12 y 12 ) e XlVl+X2V2 - e Xiy2+X2yi ^j 
Putting all together we obtain for the r.h.s. of (101): 

( ( — + xx - x 12 ~ x lXl2 yi2 ~ — + x 12 ) e Xiyi+X2y2 + (-—- Xl2 -X2 + —) e Xiy2+X2yi \ = 

x\x{ 2 y{ 2 \\y\2 2/12 J V 2/12 j/12/ J 



2^2- f(! - ^122/12) e x ^ +X2y2 - e x i y2+X2y i\ = -Mu(X,Y) (106) 
j 1/10 J 



X 122/l2 
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In order to demonstrate what the higher correlators are going to look like, we list here non-vanishing 4-point 
correlators for N = 2: 



/ 



/ UniUl) 2 [dU] = -^{(2 - 2x 12 y 12 + x 2 12 y 2 12 ) e xiyi+x ™ -2e xiy2+X2yi 

J X \ 2 V\ 2 

UuU 12 UtMi W} = ^V{(- 2 + ^122/12) e?™**™ + (2 + x l2 y l2 ) e *™+*™ 
UxxU 22 U\ x U\ 2 [dU] = -^{(2 - 2x 12 y 12 + x 2 12 y 2 12 ) e x ^ +x ^ -2e x ^ +x ^ 
U n U 22 Ul 2 U 2 \ [dU] = -J-^{(2 - x 12 y 12 ) e x ^+ x ^ - (2 + x 12 y 12 ) e x ^ +x ^ 
I UUUI) 2 [dU] = -^{2e^ + ^-(2 + 2x l2 y 12 +x 2 12 yl 2 )e^ + ^ 

J 12^12 

U 12 U 21 Ul 2 U 2 \ [dU] = -A_he x ^ +X ™ - (2 + 2x 12 y 12 + x\ 2 y 2 2 ) e x ^+ x ^ 

X \2V\2 L 



/ 
/ 



/ 



(107) 



Other non-vanishing correlators are obtained by discrete symmetries 1^2, for example, \U 22 {U\ 2 ) 2 [dU] = 
J UiiiU^) 2 [dU]. It is easy to check that the sum of the first and the second lines, which should be J UnU^ldU], 
does indeed reproduce Mn from (103). Similarly, the sum of the second and fifth line reproduce .M12, while that 
of the the second and fourth lines should - and does - vanish. Furthermore, in the limit x, y — ► expressions (107) 
turn into (5), though, as we already mentioned, this limit is not very easy to take: even in this simple situation 
(N = 2) exponents should be expanded up to the third order in their arguments. 

Actually, this last criterium considerably restricts the possible form of the 2m-point correlator: 

J u^m^^m m = 1 , (p(x 12 y 12 ) e x ™ 1+x ™ - Q(x 12 y 12 ) e *it«+*«i) (108) 

where P(t) = J2™=o ^t k and Q(t) = J2k=o Qkt k are polynomials of t = x\ 2 y\ 2 — (xi - x 2 )(y\ - y 2 ) of degree m. 
Then the condition that this expression is non-singular at t = implies that 



! 



3=0 



3 



Qk- 



(109) 



i.e. P(t) is fully defined for a given Q(t). The choice of Q(t) depends on the correlator. In particular, 
rt w rjT-n _ (~) r 



uziui r [dU] 



(£122/12)' 



ITT (E 



fe=0 



( (xi 2 yi 2 ) k ) e Xiyi+X2y2 - mle Xiy2+X2yi 



(110) 



It is clear from this example that generalization of (91) and (92) to 4-point and higher correlators should be 
straightforward, still it remains to be found. 
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